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Abstract 

Key targets for gravitational wave (GW) observatories, such as LIGO and the next 
generation interferometric detector, Advanced LIGO, include core-collapse of mas- 
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sive stars and the final stage of coalescence of compact stellar remnants. The com- 
bined GW signal from such events occurring throughout the Universe will produce 
an astrophysical GW background (AGB), one that is fundamentally different from 
the GW background by very early Universe processes. One can classify contribu- 
tions to the AGB for different classes of sources based on the strength of the GW 
emissions from the individual sources, their peak emission frequency, emission du- 
ration and their event rate density distribution. This article provides an overview 
of the detectability regimes of the AGB in the context of current and planned grav- 
itational wave observatories. We show that there are two important AGB signal 
detection regimes, which we define as 'continuous' and 'popcorn noise'. We describe 
how the 'popcorn noise' AGB regime evolves with observation time and we discuss 
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how this feature distinguishes it from the GW background produced from very early 
Universe processes. 

Key words: gravitational waves; methods: statistical; cosmology: miscellaneous 
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1 Introduction 

Astronomy has at its disposal a powerful set of technological tools to explore 
the cosmos, enabling new windows to the Universe to be opened. During the 
19th and early part of the 20th century no one envisaged the scope of astro- 
physical phenomena that the electromagnetic (EM) spectrum would reveal. 
The optical window into the Universe has been extended to include radio, 
infrared, ultraviolet, x-ray and 7-ray components of the EM spectrum, and 
neutrino detectors have been constructed and brought into successful opera- 
tion. These technology-based advances have led to many dramatic discoveries, 
including the cosmic microwave background radiation and exotic objects such 
as neutron-stars and super-massive black-holes. Gravitational wave (GW) as- 
tronomy, a new window, offers the possibility of dramatically extending our 
present understanding of the Universe. This new spectrum is presently totally 
unexplored. 

Three sets of long baseline laser-interferometer GW detectors have been, or are 
nearly, constructed. The US LIGO (Laser Interferometer Gravitational- wave 
Observatory) has started its 4th science run (S4), using two 4- km arm detec- 
tors situated at Hanford, Washington, and Livingston, Louisiana; the Han- 
ford detector also contains a 2-km interferometer. The Italian/French VIRGO 
project is commissioning a 3-km baseline instrument at Cascina, near Pisa. 
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There are detectors being developed at Hannover (the German/British GEO- 
600 project with a 600-m baseline, commissioning completed May 2006) and 
near Perth (the Australian International Gravitational Observatory, AIGO, 
initially with an 80-m baseline). A detector at Tokyo (TAMA-300, 300-m base- 
line) has been in operation since 2001. The astrophysical detection rates are 
expected to be low for the current interferometers, such as 'Initial LIGO', but 
second-generation observatories with high optical power are in the early stages 
of development; these 'Advanced' interferometers have target sensitivities that 
are predicted to provide a practical detection rate. 



The currently operating detectors are searching for both astrophysical sources 
of GWs, such as merging binary neutron-star and black-hole systems and core- 
collapse supernovae, and stochastic backgrounds composed of many individu- 



ally unresolved astrophysical sou rces at cosmological distances (jFerrari et al 



1999: 



Coward et al 



2001 



2002j). In addition to these sources, detectors may 
also be sensitive to a stochastic background of gravitational radiation from 
the early Universe, a quest that is considered by many to be the 'holy grail' 



of GW searches. 



We provide here a brief sketch of several key ideas underpinning models for 
the primordial background, but our main focus is on the other background-a 
stochastic background produced by astrophysical sources distributed through- 
out the Universe. This GW signal, as well as being of profound fundamental 
interest, provides an interesting comparison to that expected from the primor- 
dial Universe. 
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2 Cosmological background of primordial GWs 

Detection of a GW background from the early Universe would have a profound 
impact on early-Universe cosmology and on high-energy physics, opening up 
a new window to explore both the Planck epoch at 10~ 43 seconds into the 
evolution of the Universe and physics at correspondingly high energies that 
will never be accessible by any other means. Relic GWs must carry unique 
information from the primordial plasma at a graviton decoupling time when 
the Universe had a temperature ~ 10 9 K, providing a 'snapshot' of the state 
of the Universe at that epoch: cosmological GWs could probe to great depth 
in the very early Universe. 

As currently understood, mechanisms for generating GWs in the primordial 
Universe can best be described as speculative. Four popular mechanisms are 
mentioned briefly here; see Maggiore (2000, Sects. 8, 9 & 10) for a compre- 
hensive account. 

Vacuum fluctuations. In inflationary models, as the Universe cooled, it 
passed through a phase in which its energy density was dominated by vac- 
uum energy, and the scale factor increased exponentially fast. The spectrum of 
GWs that might be observed today from this period results from adiabatically- 
amplified zero-point energy fluctuations. The spectrum is expected to be ex- 
tremely broad, covering 10~ 18 — 10 8 Hz. 

Phase transitions. Phase transitions in the early Universe, particularly at 
the electro-weak scale, are potential sources of GWs. The Standard Model 
of particle physics predicts that there will be a smooth cross-over between 
states rather than a phase transition, but supersymmetric theories predict 
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first-order phase transitions, hence allowing for possible strong GW emission 
in the KT 5 - 1 Hz band. 

Cosmic strings. Grand unified theories predict topological defects; among 
these, vibrating cosmic strings, characterized by a mass per unit length, could 
potentially be strong sources of GWs. The spectrum has a narrow band feature 
peaking at 1CT 12 Hz and a broad band component in the 1CT 10 — 10 10 Hz band. 

Coherent excitations of the Universe, regarded as a 'brane' or surface 
defect in a higher- dimensional universe. The excitations could be of a 'radion' 
field that controls the size or curvature of the additional dimensions, or they 
could be of the location and shape of our Universe's brane in the higher di- 
mensions. If there was an equipartition of energy between these excitations 
and other forms of energy in the very early Universe, then the excitations 
could produce GWs strong enough for detection by advanced LIGOs. They 
could probe one or two additional dimensions of size ~ 1CT 13 - 1CT 16 m. If 
the number of extra dimensions is more than 2, much smaller scales could be 
reached (Cutler & Thorne 2002; Hogan 2000). 

2.1 Characterizing the background 

A cosmological stochastic background of GWs is expected to be isotropic, 
stationary and unpolarized. The intensity of a stochastic background of GWs 
is conventionally characterized by the dimensionless 'closure density', Qcw(f), 
defined as the fraction of energy density of GWs per logarithmic frequency 
interval normalized to the cosmological critical energy density p c c 2 required 
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to close the Universe: 



*W(/) = (i) 

p c c 2 d log / 



where pcw{f) is the energy density of the background at frequency /. In terms 
of the present value of the Hubble constant if , written as h x 100 km s" 1 
Mpc" 1 : 



Pc = « 1.9 x 10- 26 ^ kg m- 3 . (2) 



There are several other ways of describing the spectrum, including the spec- 
tral density (density in frequency space) of the ensemble average of the Fourier 
component of the metric 1 , Sh{f), with units of Hz -1 ; this allows the experi- 
mentalist to express the sensitivity of the detector in terms of a strain sensi- 
tivity, S^ 2 (f), with dimension Hz -1 / 2 . Also, the spectrum can be expressed as 
a dimensionless characteristic amplitude 2 . The relationships between h c (f), 
S h (f) and h%fLaw{f)i as derived by e.g. Maggiore (2000, Sect. 2.2), are: 



h c (f) = j2fS h (f), (3) 



K^Gw(f) = (2rr 2 /3)fX(f) = /3)fS h (f). (4) 



1 see Maggiore 2000, Sect 2.2 for a complete derivation of Sh{f ) 

2 note that the Fourier convention used here is g(f) = j^° O0 dtex.p(27rift)g(t) so 

that g(t) = J^ 0O d/exp(-27ri/t)5(/) 
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2.2 Observational bounds on Vtcw 



The integrated closure density of GWs from all frequency bands is expressed 
in dimensionless form as 



The current best limit for hlQ GW is 6 x 1(T 6 (Cutler & Thorne 2002); a value 
larger than this would have meant the Universe has expanded too rapidly 
through the era of primordial nucleosynthesis (Universe age ~ a few minutes), 
distorting the universal abundances of light elements away from their observed 
values. 

Pulsars are natural GW detectors. As a GW passes between us and the pulsar, 
the time of arrival of the pulse will fluctuate. The bound on VL GW (f) from 
pulsar timing is (Maggiore 2000, Sect. 7.3): 



A background of GWs will also cause fluctuations in the temperature of 
the cosmic microwave background. The bound from anisotropy measurements 
(Maggiore 2000, Sect. 7.2) is: 



The proposed Laser Interferometer Space Antenna (LISA) will consist of three 
spacecraft flying 5 million km apart in the form of an equilateral triangle. 
The centre of the triangular formation will be in the ecliptic plane 1 AU 
from the Sun and 20 degrees behind the Earth. The main advantage of a 
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h 2 n GW (f) < io- 8 , 



/ ~ 10" 8 Hz . 



(6) 



h 2 n GW (f) < 7 x 10- n (^ //) 2 , 1 < f/H < 30, 



(7) 
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space-based interferometer is the removal of the limitations imposed by gravity 
gradient and seismic noise, allowing a space-based detector to probe a low- 
frequency range that is not accessible to terrestrial detectors. Observations 
in the very low-frequency domain 10~ 4 Hz to 0.1 Hz are planned. Because 
Oqw(/) oc f 3 Sh(f), as in Eq. (4), an order of magnitude decrease in frequency 
results in a three order of magnitude increase in sensitivity to ^gw(/)- 

Although this increase in sensitivity is dramatic in the context of cosmological 
GW backgrounds, it introduces the problem of sensitivity to the Galactic 
population of white dwarf binaries. The signal from this Galactic source may 
prove to be a limiting factor in searches for cosmological backgrounds. 

3 Astrophysical GW backgrounds (AGB) 

3.1 Conceptual overview of the AGB 

The emission of GWs from astrophysical sources throughout the Universe can 
create a stochastic background of GWs. It could potentially be used to probe 
the Universe at redshifts z ~ 1-10 and be used as a tool to study the evolving 
star formation rate, supernova rates and the mass distribution of black-hole 
births. However, from the point of view of detecting the cosmological back- 
ground produced in the primordial Universe, the astrophysical background is 
a 'noise', which could possibly mask the relic cosmological signal. Hence, an 
understanding of AGBs is important on two fronts: first to provide fundamen- 
tal knowledge of astrophysical source evolution on a cosmological scale and 
second to differentiate this background from the early-Universe background. 
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The astrophysical sources that contribute to a GW background can be broadly 
categorized as being either continuous or burst sources. In the context of 
observing the AGB, continuous GW emission is a somewhat ambiguous ex- 
pression, as it depends on the observation time. We provide more rigorous 
definitions: 

AGB continuous sources — Astrophysical sources yielding slowly evolving 
emission with characteristic evolution time r that is very long compared with 
the observation time T. 

By choosing T 1 yr as a practical observation time, sources that could 
be considered continuous include compact binary systems in the gradual in- 
spiral phase prior to the final few minutes of their evolution and deformed 
neutron-stars. Galactic white dwarf binaries are expected to make a significant 
contribution to the AGB. 

AGB burst sources — Astrophysical sources with a GW emission time r that 
is very short compared with the observation time T. 

A typical burst source example is a core-collapse supernova, for which r could 
be of order 1CT 3 s to minutes. Others include the late in-spiral and coalescence 
stages of binary systems, neutron-star bar-mode instabilities and gamma-ray 
burst sources. 

The total event rate observed in our frame, -R to tab is obtained by integrating 
events in shells of z to the epoch when the events of the given type first started. 
For burst sources of a given type, the fraction of T that contains information 
on these events is of order i? to tai x T - Clearly, this indicates that -R to tai intro- 
duces an additional constraint for a stochastic background composed of burst 
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sources; that is -RtotaiT ~ 1. The dimensionless quantity -RtotaiT is termed duty 
cycle (DC). A DC = 1 implies a continuous signal and for a DC = 0.5, the 
signal is on average present only 50% of the time . The formal definition of 
DC that incorporates cosmological sources is: 



where r is time dilated to (1 + z)r by cosmological expansion and dR/dz, the 
variation of the event rate with z, is dependent on the cosmological model; 
the integration limit z c corresponds to the epoch when the events of the given 
type first started. 

3.2 Two AGB Components 

The AGB comprises two components: 

AGB continuous component — The contribution to the AGB from sources 
yielding a DC ^> 1, either because r is large or because i? to tai is large. The 
amplitude distribution will be Gaussian as shown by the central limit theorem, 
which states that the sum of random events from different distributions con- 
verges to a Gaussian distribution. Such backgrounds could include the gradual 
in-spiral phase of compact binary systems, e.g white-dwarf, neutron-star and 
black-hole binaries and other systems for which i? to tai is large and r is not 
extremely small. 

AGB 'popcorn noise' component — The contribution to the AGB from 
burst sources for which r is short but -Rtotai is high enough that DC ~ 1. It 
will manifest in GW data as a stochastic background comprising mostly non- 
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overlapping individual GW emissions with an amplitude distribution domi- 
nated by the spatial distribution of the sources. Possible sources include core- 
collapse supernovae, neutron-star bar-mode instabilities and GW emission as- 
sociated with gamma-ray bursts. This component is characterized by a very 
skewed distribution in GW amplitude. The distribution is found to possess 
a maximum corresponding to cosmological sources at a redshift of 2-3 and a 
long 'tail', representing more 'local', events. 

Different populations of source types (continuous and burst) will superpose to 
form a composite AGB that may be very complex. Nonetheless there could be 
single-source components that dominate at different frequencies, based on the 
simple fact that the astrophysical source types vary considerably in emission 
mechanisms. 



3.3 Spectral properties of the AGB 



The ffuence of a single source located at redshift z is given by 



1 dE GW 



where — (1 + z)r is the luminosity distance, r is the proper distance, 
dE GW /du is the angular averaged GW energy spectrum, v = (1 + z)v the 
frequency in the source frame and dR(z) is the event rate at redshift z. 

The integrated flux at the observed frequency u Q is defined by: 

F Vo = J UJR(z) (10) 
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Using the definitions above, a stochastic gravitational wave background of 
astrophysical origin, eq.(l) can be expressed as 

tt GW {vo) = —v F Vo (11) 

C Pc 

Different AGB models have been investigated (see Regimbau 2005 for a com- 
plete review) by using several astrophysical sources of GWs. On the one hand, 
distorted black-holes (Ferrari et al. 1999; de Araujo et al. 2000) are examples 
of sources able to generate a shot noise (time interval between events large in 
comparison with duration of a single event), while supernovae or hypernovae 
(Blair et al. 1996; Coward et al. 2001; Coward et al. 2002 ; Howell et al. 2004; 
Buonanno et al. 2004) are expected to produce an intermediate 'popcorn' noise 
i.e. between a shot noise and continuous signal. On the other hand, the con- 
tribution of tri-axial rotating neutron-stars (Regimbau & de Freitas Pacheco 
2001), including magnetars (Regimbau & de Freitas Pacheco 2005a), consti- 
tutes a truly continuous background. Figure (1), from Reigimbau 2005, shows 
different types of AGB spectra along with spectra of a primordial cosmological 
background. 

3.4 Continuous AGB to single event 

The AGBs composed of burst sources discussed above are assumed to consist of 
a continuous (or nearly continuous) background of GWs that would be orders 
of magnitude below the noise level of any single Advanced LIGO-type detector. 
Qualitatively one expects that the AGB will be dominated by sources typically 
3-4 orders of magnitude fainter in luminosity than sources in relatively nearby 
galaxies at 10-100 Mpc (z « 0.002-0.02). 
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There is a crucial difference between a primordial cosmological GW back- 
ground and the AGB: the primordial background is long-term stationary over 
any practical observation time T but, in contrast, the AGB signal over T will 
fluctuate according to dR/dz. To explain the latter point, consider GW events 
associated with stellar core collapses. Assume these events occur at a local rate 
of 1 yr^ 1 out to 15 Mpc and produce a signal above the detectable threshold 
in an idealized GW detector at this source distance. For T <C 1 yr the detector 
will be dominated by a roughly stationary background from sources at z ~ 
1-2, well below the detectable threshold. But, as T — > 1 yr, it will be highly 
probable that the GW detector will record a local event. What is observed as 
a local event can be interpreted as the low probability tail of dR/dz or the 
local component of the AGB. 

The low probability tail contains likely sources that may be detectable as sin- 
gle events by LIGO and Advanced LIGO over observation times of about 1 yr: 



A significant part of current detection strategies are focussed on trying to de- 
tect the low probability tail of the popcorn noise distribution. 



The high- ,2 sources will accumulate as unresolved signals in the data and with 
no prior knowledge of this background, it will manifest as an excess noise in 
the detector. The optimal signal detection method is cross-correlating two or 
more detectors. 
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4 The probability event horizon 



For any GW detector and GW source type, one can define a 'detectability 
horizon' centred on the detector and surrounding the volume in which such 
events are potentially detectable; the horizon distance is determined by the 
flux limit of the detector and the source flux. A second horizon, defined by 
the minimum distance for at least one event to occur over some observation 
time, with probability above some selected threshold, can also be defined. 
We call this the probability event horizon, or the PEH 3 . It describes how an 
observer and all potentially detectable events of a particular type are related 
via a probability event distribution encompassing all such events. The PEH 
is a way of quantitatively describing how the popcorn noise component of the 
AGB will manifest in GW detector data. 

The idea can be illustrated using a static Euclidean universe, assuming an 
isotropic and homogenous distribution of events. For a constant event rate 
r per unit volume, the mean cumulative event rate in a volume of radius 
r is R(r) = (4/3)7rr 3 r . The events are independent of each other, so their 
distribution is a Poisson process in time: the probability for at least one event 
to occur within the radius r is: 



e R being the probability of zero events occurring. For a chosen value of p, 
one solves the equation p = e for r, the corresponding radial distance. The 
3 see Coward k, Burman (2005) for a more detailed description of the PEH in a 
broader astronomical context. The authors apply it to gamma-ray burst red-shift 
data to demonstrate how it can be used to constrain the local rate density of GRBs. 



p(n > l;R{r),T) = 1 -e 



-R{r)T 
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PEH is defined as: 

r PEH (T) = (3iVe /4 7rro )l/3 T -l/3 



(13) 



The speed at which the horizon approaches the observer — the PEH velocity- 
is obtained by differentiating r e (T) with respect to T: 



PEH 



(T) = (N e /36nr ) 1/3 T~ 4/3 . (14) 



The PEH concept is a useful tool for modelling the popcorn component of 
the AGB, such as the background from double neutron-star (DNS) mergers. 
A PEH model can in principle be used to constrain the local rate density of 
events in GW detector data, as shown by Howell et al. (2006). It also provides 
a way of describing the detectability of the AGB for a certain GW detector 
sensitivity, such as LIGO or Advanced LIGO (Coward et al. 2005). Figure 2 
plots the PEH for DNS mergers, assuming the events are roughly proportional 
to the evolving star-formation rate (SFR). Coward et al. (2005) justify this 
assumption by showing that a 1-2 Gyr merger time does not significantly alter 
the distribution of mergers in redshift, especially for small z. The plot shows 
how the horizon approaches the detectability volume encompassing the LIGO 
and Advanced LIGO detectors. The time evolution of the AGB described by 
the PEH shows how the rate density of events goes from a high rate of distant 
low-luminosity events to rare higher luminosity events. 



5 Detection regimes of the AGB 

If the number of GW sources is large enough that the time interval between 
events, as observed in our frame, is small compared to the duration of a single 
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event (DC » 1), the combined GW emissions will produce a continuous 
background. The central limit theorem predicts that this regime of the AGB 
will yield a Gaussian signal. For this case the optimal detection strategy is to 
cross-correlate the output of two (or more) detectors, assuming they have in- 
dependent spectral noises. Matched filtering may also be used here to optimise 
the cross-correlation. 

In the AGB regime where the source rate density is small and the time interval 
between events is long compared to the duration of a single event (DC « 1), 
the sources are resolved and the nearest may be detected by 'burst' data 
analysis techniques. The optimal method is matched filtering but this filter 
can only be optimal if the exact shape of the signal is known. Robust generic 
methods based on coincident analysis are also being developed to search for 
events where the exact form of the signal is not known (see for instance Arnaud 
et al. 1999; Pradier et al. 2001; Sylvestre 2002; Chassande Mottin & Pai 2005). 

In the intermediate 'popcorn noise' regime, where the time interval between 
events is of the same order as the duration of a single event, the waveforms 
may overlap but the amplitude distribution will not be Gaussian. The am- 
plitude distribution of this background signal will contain information on the 
spatial distribution of the sources. This point is highlighted by using the PEH 
description for 'popcorn noise' : the amplitude of the individual events scales 
with distance and the distance to the largest-amplitude event is reducing as 
a function of observation time, as shown in figure 2. The time dependence 
of the signal is important and data analysis techniques in the frequency do- 
main, such as the cross correlation statistic, do not utilze this feature of the 
signal. New data analysis techniques are currently under investigation, such 
as the search for anisotropies (Allen and Ottewill, 1997) which can be used 
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to create a map of the GW background (Cornish 2001; Ballmer 2005), the 
maximum likelihood statistic (Drasco & Flannagan 2003), or by utilizing the 
PEH concept described in the previous section (Howell et al. 2006). 

The different regimes of the AGB are being studied using population synthesis. 
Regimbau & de Freitas Pacheco (2005b) performed numerical simulations of 
the cosmological population of coalescing neutron-star-neutron-star binaries. 
Figure 3 shows that events located beyond the critical redshift = 0.23 
produced a continuous background, while those in the redshift interval 0.027 < 
z < 0.23 produce a nearly continuous (popcorn noise) signal. The signal is 
discrete, i.e. shot noise like, for events occurring closer than z = 0.027. The 
GW closure density, Q gw , for the continuous component, has a maximum at 
670 Hz with an amplitude of 1.1 x 10~ 9 , while the popcorn noise component 
has an amplitude about one order of magnitude higher with maximum at 1.2 
kHz. That study highlights the importance of the popcorn noise regime to the 
AGB. 

6 Discussion 

We show that the AGB can be classified into several detection regimes. For 
continuous GW backgrounds, cross correlating multiple detectors is the opti- 
mal search method, assuming that the noise in each detector is not correlated. 
Examples of GW source populations where this technique is appropriate in- 
clude both Galactic binaries and massive cosmological binary systems. The 
central limit theorem predicts that the sum of the GW signal from these 
sources produces an amplitude distribution that is Gaussian. Because of the 
inverse distance dependence for the GW amplitude, this continuous regime 
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may be dominated by a small number of relatively close systems, implying 
that matched filtering will be the optimal search method. This is an example 
of how the low probability tail of the GW event rate density corresponding 
to nearby sources, dominates the signal over 'long' observation times. 

We identify another AGB regime, where the GW signal composed of 'burst 
souces' contains information on the source-rate distribution. An example of the 
sources contributing to such a background include the GW emission associated 
with gamma-ray bursts and massive stellar collapse. Even if the cosmological 
event rate for these sources is high enough to produce a truly continuous GW 
background, there is a horizon distance where the signal will be dominated 
by nearby events. Hence it is apparent that there is a transition from a small 
amplitude continuous unresolved signal to higher-amplitude discrete events. 

The detectability of the AGB is determined by the detector sensitivity which 
sets the detection horizon. If the detector horizon extends to distances cor- 
responding to event rates for which the received signal is nearly continuous, 
i.e. to the popcorn noise regime, it is possible to probe the source distribution 
and its association with the evolving star formation rate. Because the present 
detector horizon for LIGO is of some tens to hundreds of Mpc for the strongest 
proposed GW sources, the AGB regime currently probed is dominated by rare 
single events corresponding to relatively long observation times, as shown by 
the PEH curves in Figure 2. As the sensitivity of LIGO improves and the next 
generation of GW observatories come on line, it is plausible that the AGB 
can be explored for the first time, providing new insight into both the sources 
generating the AGB and the evolutionary history of the Universe. 
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Fig. 1. The dimensionless gravitational wave closure density, Fl gw , plotted as a 
function of observed frequency for potentially significant contributions from astro- 
physical GW sources: continuous background (continuous line) such as from neu- 
tron-star-neutron-star coalescence (Regimbau &; de Freitas Pacheco 2005b), rotat- 
ing pulsars or magnetars (Regimbau & de Freitas Pacheco 2005a) and shot noise 
(dotted line) such as black-hole ringdown and bar-mode instabilities in young, fast 
rotating neutron-stars (Regimbau 2005b). Theoretical limits for several primordial 
GW background models are plotted (horizontal lines) and the detection sensitivity 
from correlating two Advanced LIGO type interferometers is shown (dashed line). 
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Fig. 2. The PEH evolution for double neutron-star mergers, from Coward et al. 
(2005), as a function of observation time. Galactic merger rate limits are very un- 
certain, but we use 292 Myr" 1 (dashed curve) and 17 Myr -1 (solid line) as examples, 
and we assume the merger rate evolution following the evolving SFR. The horizon- 
tal lines show the sensitivity limits for Advanced and Initial LIGO and the PEH 
curves assume an idealized optimal detection filter. 
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Fig. 3. Same as Figure 1, but showing the different regimes of the AGB for double 
neutron-star mergers occurring beyond = 0.23 (bold continuous curve). The 
regime with duty cycle DC > 10, corresponds to sources beyond z = 1.05 (grey 
continuous curve). The 'popcorn' noise regime, with DC = 1, corresponds to sources 
in the redshift interval 0.027 < z < 0.23 (solid line), and the 'shot' noise regime, 
with DC < 0.1, corresponding to sources at z < 0.027. The GW closure density is 
dominated by the popcorn and shot noise contributions. 
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